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<N ; ABSTRACT 

o . 

O ! We present results of long-slit spectroscopy in several positions of the Orion 

0> I nebula. Our goal is to study the spatial distribution of a large number of neb- 

I— i! ular quantities, including line fluxes, physical conditions and ionic abundances 

p ^. at a spatial resolution of about 1". In particular, we have compared the O"'""'" 

Q ! abundance determined from collisionally excited and recombination lines in 671 

^ ! individual ID spectra covering different morphological zones of the nebula. We 

! find that protoplanetary disks (proplyds) show prominent spikes of Te([N II]) 

! probably produced by collisional deexcitation due to the high electron densities 

^ i found in these objects. Herbig-Haro objects show also relatively high Te([N II]) 

! but probably produced by local heating due to shocks. We also find that the 

^ I spatial distribution of pure recombination O II and [O III] lines is fairly similar, 

^ J I in contrast to that observed in planetary nebulae. The abundance discrepancy 

O I factor (ADF) of 0++ remains rather constant along the slit positions, except in 

I some particular small areas of the nebula where this quantity reaches somewhat 

^ I higher values, in particular at the location of the most conspicuous Herbig-Haro 

I objects: HH 202, HH 203, and HH 204. There is also an apparent slight increase 

^ I of the ADF in the inner 40" around 6^ Ori C. We find a negative radial gradient 

5^ I of Te([0 III]) and Te([N II]) in the nebula based on the projected distance from 

6^ Ori C. We explore the behavior of the ADF of O"*"^ with respect to other 
nebular quantities, finding that it seems to increase very slightly with the elec- 
tron temperature. Finally, we estimate the value of the mean-square electron 
temperature fluctuation, the so-called parameter. Our results indicate that 
the hypothetical thermal inhomogeneities -if they exist- should be smaller than 
our spatial resolution element. 

Subject headings: ISM: abundances -H II regions- ISM: individual: Orion nebula 
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Introduction 



The analysis of tlie spectrum of H II regions allows to determine the chemical compo- 
sition of the ionized gas phase of the interstellar medium from the Solar Neighbourhood to 
the high-redshift Universe. Therefore, it stands as an essential tool for our knowledge of the 
chemical evolution of the Universe. In photoionized nebulae, the abundance of the elements 
heavier than He is usually determined from collisional excitation lines (hereinafter CELs), 
whose intensity depends exponentially on the electron temperature, Te, of the gas. It was 
about 20 years ago when the first determinations of C^^ abundance from the intensity of 
the weak recombination line (hereinafter RL) of C II 4267 A were available for planetary 
nebulae (PNe). The comparison of the abundance obtained from C II 4267 A and from the 
CELs of this ion in the ultraviolet (UV) showed a difference that could be as large as a orde r 



of magnitude in some objects (e.g. lFrenchlll983l : iRola &: Stasihskalll994j : iMathis fc Liijll999l ) . 
Peimbert et al.l (119931 ) were the first in determinig the 0"'"+ abundance from the very weak 



RLs, obtaining the same qualitative result: the abundances obtained from RLs are higher 
than those determined making use of CELs. This observational fact is currently known as 
the "abundance discrepancy" (hereinafter AD) problem. In the last years, our group has 
obtained a large dataset of intermediate and high resolution spectrosc opy of Galactic and 



extragalactic H II regions using medium and larg e aperture telescopes (lEste 



ban et al. 2002 



2OO5I : ICarcfa-Roias et all hooi boosl . bood . I2OO7I : iLopez-Sanchez et al.ll2007h . The general 
result of these works is that the 0"'""'"/H"'" ratio calculated from RLs is between 0.10 and 0.35 
dex higher than the value obtained from CELs in the same objects. The value of the AD that 
we usually find in H II regions is rather similar for all objects and ions and is much lower than 
the most extreme values found in PNe. The results for H II regions obtained by our group 
are fairly different that those found for PNe, and seem t o be consis t ent w ith the predictions 
of the temperature fluctu ations paradigm formulat e d by PeimbertI (119671 ). as it is argued in 



Garcia- Roj ad (120061 ) and iGarcfa-Rojas fc EstebanI (120071 ) . In the presence of temperature 



fluctuations (parametrized by the mean square of the spatial variations of temperature, the 
so-called t"^ parameter) the AD can be naturally explained because the different temperature 
dependence of the intensity of RLs and CELs. The existence and the origin of temperature 
fluctuations are sti ll controversial problems an d a c hallenge for our unders tanding of ionized 
nebulae. Recently. iTsamis &: PequignotI (120051 ) and lStasihska et al.l (120071 ) have proposed an 
hypothesis to the origin of the AD, which is based on the presence of cold high-metallicity 
clumps of supernova ejecta still not mixed with the ambient gas of the H II regions. This 



^Based on observations made with the 4.2m WiUiam Herschel Telescope (WHT) operated on the island 
of La Palma by the Isaac Newton Group in the Spanish Observatorio del Roque de los Muchachos of the 
Instituto de Astrofisica de Canarias. 
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cold gas would produce most of the emission of the RLs whereas the ambient gas of normal 
abundances would emit most of the intensity of CELs. 

Our group is interested in exploring on what variable or physical process the AD de- 
pends from different approaches. One of the most promising is based on the study of the 
behaviour of this magnitude at small spatial scales, something that has still not been ex- 
plored in depth in nearby bright Galactic H II regions. In this paper, we make use of deep 
intermediate-resolution long-slit spectroscopy of the Orion nebula to study the dependence of 
the AD with respect to different nebular parameters: electron temperature and density, local 
ionization state of the gas, presence of high velocity material, and its correlation with differ- 
ent morphological structures (proplyds, ionization fronts, globules, Herbig-Haro objects), in 
H II regions. 

The spatial distribution of the physical condit ions in the Orion nebula has been in- 
vestigated by several authors. [Baldwin et al.l (Il99ll ) obtained the density and temperature 
distribution in 21 and 14 points, respectively, along a 5' line west of 6^ Or i C, finding a density 
gradient that decreases to the outskirts of the nebula and a constant Tg. IWalter et al.l ( 1l992l ) 
determined electron densities and temperatures and chemical abundances for 22 regions of 
the Orion nebula. Using also data from the literature, these authors find radial gradients of 
the physical condit ions, but wit h a pos itive slope in the case of the temperature determined 
from [O III] lines. iPogge et al.l (119921 ) obtained Fabry-Perot images of the inner 6' of the 
nebula covering several bright CELs and taken with an average seeing of about I'.'S. Those 
authors present a density map obtained from the ratio of the [S II] doublet confirming the 
presence of a density gradient that reaches its highest point immediately south-southwest of 
the Trapezium stars, and some locali zed density en h ancem ents in the Orion bar and some 
Herbig-Haro objects. Very recently, ISanchez et al.l (120071 ) have obtained an integral field 
spectroscopy mosaic of an area of 5' x 6' of the center of the Orion nebula, with a spatial 
resolution of 2'! 7. T he electron density map they obtain is consistent with that obtained by 
Pogge et al.l (Il992l) but richer in sub structures, some of them possibly associated to Herbig- 
Haro objects. ISanchez et al.l (120071 ) also obtain a electron temperature map (derived from 
the line ratio of [N II] lines) that shows clear spatial variations, which rise near the Trapezium 
and drop to the outer zones of the nebula. However, an important drawback of the tempera- 
ture map of Sanchez et al. is that is based on non-flux calibrated spectra and p ossible effects 



due to variations in the dust extinction distribution cannot be disregarded. lO'Dell et al. 



( 120031 ) obtained a high spatial resolution map of the electron temperature -derived from the 
line ratio of [O III] lines- of a 160" x 160" field centered at the southwest of the Trapezium. 
The data were obtained from narrow-band images taken with the WF PC of the HST. Al- 
though they do not find a substantial radial gradient of Te in the nebula, lO'Dell et al.l (120031 ) 
report the existence of small-scale temperature variations down to a few arcseconds com- 
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patible with the v alues of the temperat u re fluctuations para meter calculated from the AD 
determinations by lEsteban et al.l (120041 ) . iRubin et al.l (120031 ) obtained HST /STIS long-slit 
spectroscopy at several slit positions on the Orion nebula analysing the electron tempera- 
ture and density spatial profiles with resolution elements of O'.'S x O'.'S. These last authors 
do not find large-scale gradients of the physical conditions along the slits but a relatively 
large point-to-point variation and some correlation of such variations with several small-scale 
structures. 

The spatial mapping of t he A D factor has been perfor r ned in few i onized nebulae but 



l argely for PNe. iLiu et all (l2000f ). iGarnett fc DinersteinI (120011 ). and iKrabbe fc Copetti 
(120061 ) have found significant differences in the spatial profiles of the 0+"'"/H+ ratio de- 



rived making use of RLs and CELs suggesting the presence of chemic al inhomogene i ties or 
additional mechanisms for producing the O II lines in these objects. iTsamis et al.l (120031 ) 
have performed the only available study so far of the spatial distribution of the AD factor 
in an H II region: 30 Doradus. However, considering the extragalactic nature of this object 
and the spatial sampling of 3^^5 used by those au thors, their final spatial resolution is very 
low -about Ipc. In any case, iTsamis et al.l (l2003l ) find a rather constant AD factor along the 
zone covered with their observations, a quite different behavior than that observed in PNe. 

In §§ |2]and|3]of this paper we describe the observations, the data reduction procedure 
and the aperture extraction and measurement of the emission lines. In § H] we derive the 
physical conditions and the ionic abundances from both kinds of lines: CELs and RLs. In 
§ we present and discuss the spatial profiles of the physical conditions, hne fluxes, and 
the abundance discrepancy factor along the slit positions. In § |6] we discuss the large-scale 
radial distribution of some nebular properties along the nebula. In § |7] we explore possible 
correlations between the AD and different nebular parameters. In §|8]we address and estimate 
the possible temperature fluctuations inside the nebula. Finally, in § IH] we summarize our 
main conclusions. 



2. Observations, Data Reduction, and Extraction of ID Spectra 

Intermediate-resolution spectroscopy was obtained in 2002 December 27 with the ISIS 
spectrograph at the 4.2m William Herschel Telescope (WHT) in Observatorio del Roque de 
los Muchachos (La Palma, Spain). Two different CCDs were used at the blue and red arms 
of the spectrograph: an EEV CCD with a conflguration 4096 x 2048 pixels with a size of 
13.5 /xm per pixel in the blue arm and a Marconi CCD with 4700 x 2148 pixels with a pixel 
size of 13.5 fim in the red arm. The dichroic prism used to separate the blue and red beams 
was set at 5400 A. The slit was 3'.7 long and 1'.'03 wide. Two gratings were used, the R1200B 
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in the blue arm and the R316R in the red arm. These gratings give reciprocal dispersions 
of 17 and 62 A mm~^, and effective spectral resolutions of 0.86 and 3.81 A for the blue and 
red arms, respectively. The blue spectra cover from AA4198 to 5048 A and the red ones from 
AA5370 to 8690 A. The spatial scale is 0'.'20 pixel"^ in both arms. The average seeing during 
the observations was ~1'.'2. 

We observed 5 slit positions covering different zones of the nebula and different position 
angles (see Figure [1]). These positions were chosen in order to cover different morphological 
structures as proplyds (158-323, 158-326, 159-350, 170-337, and 177-341), Herbig-Haro ob- 
jects (HH 202, HH 203, HH 204, HH 529, and HH 530) and the Orion bar. Due to the high 
surface brightness of the nebula, a large number of individual short exposures were taken in 
each slit position and spectral range in order to achieve a good signal-to-noise ratio in the 
faint C II and O II RLs and to avoid saturation of the brightest emission lines. The journal 
of observations can be found in Tabled! Note that the 5 slit positions are numbered as 1, 3, 
4, 5, and 6, actually there was not a slit position number 2. 

The spectra were wavelength calibrated with a CuNe-|-CuAr lamp. The correction for 
atmospheric extinction was performed using the average curve for continuous atmospheric ex- 
tinction at Roque de los Muchachos Observatory. The absolute flux calibration was achieved 
by observations of the standard stars Feige 15, Feige 110, H600, and Hz 44. All the CCD 
frames were reduced using the standard TWODSPEC reduction package to perform 

bias correction, flat-fielding, cosmic-ray rejection, wavelength and flux calibration, and sky 
subtraction. 

The extraction of ID spectra was done automatically through an I RAF script, using 
the apall task. First, we traced the apertures interactivelly by selecting the brightest object 
(star, proplyd or Herbig-Haro object) in each 2D spectrum. In all cases we adjusted a third- 
order spline function and obtained a typical rms of the fit between 0.05 and 0.1 pixels. 
In the following step, we defined the apertures to extract for each 2D spectra by using 
the coefficients obtained in the aperture tracing. For all the slit positions, we extracted 
apertures of 6 pixels size in the spatial direction, that corresponds to an angular scale of 1'.'2 
-the average seeing during the night. Therefore, each aperture covers a size of 1'.'2 x 1'.'03 in 
the spatial a nd spectral directions respectively. At the distance of the Orion nebula (450 pc. 



O'Delll 120011 ). 1" corresponds to a linear size of 0.0022 pc, 6.8 x 10"^^ cm, or 450 A.U. The 
slit center in the red arm is some pixels displaced with respect to the slit center in the blue 
arm; this effect has been corrected in the extraction procedure, ensuring the same spatial 



■^IRAF is distributed by National Optical Astronomical Observatories, operated by the Associated Uni- 
versities for Research in Astronomy,Inc., under contract to the National Science Foundation 




Fig. 1. — Our slit positions over a mosaic of a cornbinat ion of WFPC2 images of the Orion 
nebula taken with different filters (jO'Dell fc Wonglll996l ). 
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covering in the blue and red ranges. We have discarded the apertures located near the edges 
of the CCD, resulting in a final number of 154 apertures extracted -individual ID spectra- 
for each slit position except numbers 1 and 4. In these two cases a star -6^ Ori A- fell into 
the slit, and we discarded the apertures contaminated by stellar emission (11 apertures in 
position 1 and 12 in position 4). In addition, the last 17 apertures at the northwest edge of 
position 4 were also discarded because the temperature sensitive [O III] 4363 A fine was not 
detected due to the faintness of the spectra. The final total number of apertures extracted 
was 730. 

For each slit position, we extracted additional ID spectra collapsing the whole slit -the 
whole extension of the 154 individual apertures- but excluding different particular zones: a) 
only those apertures contaminated by stellar emission -these will be designated as "whole 
slit" spectra, or b) the same zones as in a) but also those apertures covering proplyds, 
Herbig-Haro objects, or having a very low surface brightness -these spectra are designated 
as "background gas" . In Table [2], for each slit position, we summarize the number of apertures 
excluded due to stellar emission contamination or non detection of the [O III] 4363 A line 
(upper row), the total number of usable apertures (middle row), and the apertures that we 
consider representative of the "background gas" (lower row). 

3. Emission Line Selection, Flux Measurements, and Reddening Correction 

The emission lines considered in our analysis are indicated in Table O These lines were 
selected according to the following criteria: 

• HI lines -Ha, H/3 and H7-, which are used to compute the reddening correction and 
to re-scale the line intensity ratios of the red spectral range with respect to the blue 
one. 

• Ratios of CELs of various species, which are used to compute physical conditions -such 
as the auroral lines [O III] 4363 A or [N II] 5754 A to derive Tg or [S II] lines to derive 
rig- and ionic abundances. 

• Faint recombination lines of C II and O II, which are used to derive the C"*""*" and 
0^+ abundances and to compute the abundance discrepancy factor for 0++ (via a 
comparison with the O^"^ abundances from CELs). 

• Some lines which are blended with other lines of interest. 
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Table 1. Journal of observations 













Spectral 


Exposure 


Slit 






P.A. 


AA 


Resolution 


Time 


Position 


R.A.'' 


Decl.'' 


(dcg) 


(A) 


(A pix-l) 


(s) 


1 


05 35 15.0 


-05 23 04 





4198-5048 


0.23 


30x60 










5370-8690 


0.84 


40x30 


3 


05 35 22.4 


-05 25 09 


147 


4198-5048 


0.23 


15x150 










5370-8690 


0.84 


25x60 


4 


05 35 12.4 


-05 22 44 


107 


4198-5048 


0.23 


31x60 










5370-8690 


0.84 


40x30 


5 


05 35 15.8 


-05 23 33 


50 


4198-5048 


0.23 


20x100 










5370-8690 


0.84 


39x30 


6 


05 35 15.3 


-05 23 38 


72 


4198-5048 


0.23 


20x100 










5370-8690 


0.84 


41x30 



"■Coordinates of the slit center 



Table 2. Number of apertures extracted 





Pos 1 


Pos 3 


Pos 4 


Pos 5 


Pos 6 


No. of excluded apertures*' 


11 




29 






No. of extracted apertures'^ 


143 


154 


125 


154 


154 


"Background gas" zones'^ 


1-60,72-84,90-143 


85-154 


25-50,75-85 


1-57,62-74,79,154 


1-40,45-63,70-154 



''Contaminated by stellar emission or non detection of the [O III] 4363 A line. 
''Summed for producing the "whole slit" spectra. 
'^Summed for producing the "background gas" spectra. 
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Line fluxes were measured applying a single or a multiple -in the case of line blending- 
Gaussian profile fit procedure except in some apertures of positions 3 and 4 where the 
complex velocity field of the Herbig-Haro objects affects the line profiles and the Gaussian 
fit was not feasible. In these cases, the line intensities were measured by integrating all 
the flux included in the line profile between two given limits and over a local continuum 
estimated by eye. All the measurements were made with the SPLOT routine of the IRAF 
package and using our own scripts to automatize the process. 

To accurately compute the line fluxes we need to define the adjacent continuum of each 
line. This was done also by using the SPLOT routine. For each selected line, we define two 
small spectral zones -one at each side of the line- and located as close as possible of the line 
and free of any spectral feature. SPLOT fits the continuum between both zones and obtains 
the flux of the line. This automatic procedure was tested in randomly selected lines where we 
manually measured the flux using a local continuum estimated by eye. In general, there is a 
good agreement between both kinds of measurements within the adopted uncertainties. The 
observational errors ass ociated with the line flux measurements was determined following 



Castellanos et al.l (120021 ) from the expression (T\ = ac 

ArV2[i + EW/(A^A)]^/^ where 0"! is 
the error in the line flux, Uc represents the standard deviation of the continuum close to 
the measured emission line, N is the number of pixels used in the measurement of the line 
flux, EW is the line equivalent width, and A is the wavelength dispersion in A pixel"^. 
The final uncertainty of the line intensity ratios are estimated to be tipically: about 1% 
if the ratio F(A)/F(H/3)>0.1, about 2% if 0.01< F(A)/F(H/3) <0.1, about 10% if 0.005< 
F(A)/F(H/5)<0.01, about 25% if 0.001< F(A)/F(H/3)<0.005 and about 40% if 0.0001< 
F{X)/F{Rf3)< 0.001. We do not consider lines weaker than 0.0001 x F(H/?) (see below). 

All the line intensities of a given aperture have been normalized to a particular H I 
recombination line present in each wavelength interval. For the blue spectra, the reference 
line was H/3, and for the red one, the reference was Ha. To produce a final homogeneous 
set of line intensity ratios, all of them were re-scaled to H/5. The re-scaling factor used in 
the red spectra was the theoretical Ha;/II/5 ratio for the physical conditions of Tg = 10000 
K and ng = 1000 cm~^. 

As it can be expected, weak lines were not detected in all the apertures. To avoid bad 
weak line measurements, we imposed three criteria to discriminate between real features and 
noise. The criteria are the following: 



• Line intensity peak over 2.5 times the sigma of the continuum {2.5a): The more usual 
criterion of 3a was not adopted because several important lines clearly detected were 
discarded using that constraint. 
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• FWHM(A) > 1.5 X FWHM(H I) or FWHM(A) < FWHM(H I)/1.5: These inequalities 
were used to discriminate between true emission lines and spurious features. 

• /(A) < 0.0001 x/(H/5). This intensity was near the detection limit of our observations. 



The reddening coefficient, c (H/3), was obtained by fit ting the observed H7/H/5 ratio to 
the theoretical one predicted by Storev fc Hummerl ( 19951) for the neb ular conditions deter- 
mined in the slit position observed by lEsteban et al. kooi \. Following [j^teban et all Jl998h 



we ha ve used the reddening function, /(A), normalized at H/3 derived by lCostero fc Peimbert 



(I1970I ) for the Orion nebula. In Figure [21 as an example, we show the spatial distribuition of 
c(H/3) obtained for the apertures of the slit position 6. The zone of the largest reddening at 
the left half of Figure [2] correspond s to th e nor theast portion of slit position 6. The redden- 
ing maps of lO'Dell fc Yusef-Zadehl (I2OOOI ) and ISanchez et al.l (120071 ) also show larger values 
at this zone, which is in the vicinity of the so-called Dark Bay. The typical uncertainty of 
c(H/5) is estimated to be about 0.05. 



4. Physical Conditions and Chemical Abundances 
4.1. Electron Temperatures and Densities 



Nebular electron temperatures, T^, and densities, n^, have been derived from the usua. 1 



CEL ratios, using the IRAF t ask temden of the packa ge nebular ( jShaw &: Dufourl Il995l ) 
with updated atomic data (see iGarcia-Rojas et al.ll2005l ). We have computed Ug from the 
[S II] 6717/6731 line ratio and from the nebular to auroral [O III] (4959-1-5007) /4363 
and [N II] (6548-1-6584) /5754 line ratios. The spatial distributions of the physical conditions 
are presented and discussed in §5.11 Although we include the [CI III] doublet in our set of 
selected lines we do not use the lines in the analysis because the number of apertures with 
good determinations of the density sensitive line ratio is rather low. 

The methodology for the determination of the physical conditions was the following: 
an initial of 10000 K was assumed in order to derive a first approximation to ne([S II]); 
then, the obtained rie was used to compute Te([0 III]) and Te([N II]), and finally we iterated 
until convergence to compute the finally adopted values of rig and Tg. Uncertainties in the 
physical conditi ons were computed by propagating the errors i n the analytical expression o f 
rie computed by lCastaneda et al.l (Il992l ) and that of Tg given bv lOsterbrock &: Ferlandl (|2006l) 
(their equations 5.4 an 5.5). Although the expression derived by ICastaneda et aL (1992) is 
only valid to a limited range of densities lower than 10^ cm~^, it seems adequate for simply 
estimating the error propagation due to uncertainties in the computed temperatures and line 
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Fig. 2. — Spatial distribution of the reddening coefficient, c(H/3), along slit position 6. 
Positional measurement along the slit goes from northeast to southwest (see Figured]). The 
dashed horizontal line represents the value of c(H/5) obtained from the integrated spectrum 
along the whole slit. The typical error bar is included. 
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ratios. 



4.2. Ionic Abundances from CELs and RLs 

Ionic abundances of N^, O^, O^^, S^, S^^, and Ar^+ have been derived from CELs 
making use of the IRAF task ionic of the package nebular. We have assumed a two-zone 
scheme and = 0, adopting Te([N II]) for ions with low ionization potential (N"*", 0+ and 
S^) and Te([0 III]) for ions with high ionization potential (0++, and Ar++). The errors 
in the ionic abundances are the quadratic sum of the independent contributions of Ue, Te, 
and line flux uncertainties. 

On the other hand, the high signal-to-noise of the spectra has permitted us to detect 
and measure pure RLs of O II and C II in most of the apertures (see Figure [3]). These lines 
have the advantage that their relative intensity with respect to H I lines depend weakly on 
Te and ne, avoiding the problem of the presence of temperature variations along the line of 
sight, that can actually affect the abundance determinations from CELs, which emissivities 
depend strongly on the Tg of the nebula. 

Let /(A) be the intensity of a RL of a element X, i times ionized, at wavelength A; then 
the abundance of the ionization state i + 1 of element X is given by 

A(X-+i) _ A(A)«e//(H/?) /(A) 

A(H+) 4861 ae//(A) /(H/?) ' ^ ^ 

where ae//(A) and ae//(H/3) are the effective recombination coefficients for the line and H/3, 
respectively. The Q;e//(H/3) /«£//( A) ratio is almost independent of the adopted temperatures 
and densities. 



Following lEsteban et al.l (119981 ) we have considered the abundances obtained from the 
intensity of each individual line of multiplet 1 of O II and the abundances from the estimated 
total intensity of the multiplet. This last quantity is obtained by multiplying the sum of the 
intensities of the individual lines observed by the multiplet correction factor, defined as: 

_ _ J2all ^ij /r,\ 
^cf - ^ , (2) 

where Sij is the theorical line strength, which are constructed assuming that they are propor- 
tional to the population of their parent levels assuming LTE computations predictions. The 
upper sum runs over all the lines of the multiplet and the lower sum runs over the observed 
lines of the multiplet. 
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Fig. 3. — Portion of the blue spectrum of an aperture extracted from slit position 6 and 
extending from the positional measurements 151.4 to 152.6 arcseconds (see Figure E]). This 
is a representative example of the average quality of our ID spectra. 
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The O^^ and C^'^ abundances from RLs have been calculated using the representative 
Tg of thes e ions -Tp([0 III])- and the effective recombina tion coefficients available in the 
literature (IStoreylll994l for O II assuming LS coupling, and lDavey. Storey. &: Kisieliuj 12000 
for C II). The 0"*"+ abundance has only been computed when at least three lines of multiplet 
1 were measured in a given ID spectrum. The final number of apertures with determinations 
of the 0^"'"/H+ ratio obtained from RLs was 671, a 92% of the total number of available ID 
spectra. NLTE cor r ections are not taken into account for d eriving the O"*"^ abundances (see 



Tsamis et al.l l2003l : iRuiz et al.l l2003l : iPeimbert et al.l |2005| ) considering that we use several 



lines of the multiplet and the high densities -between 4000 and 6000 cm ^ (e.g. lEsteban et al. 



19981 . |200J)- of the Orion nebula. 



5. Spatial Profiles along the Slit Positions 

5.1. Physical Conditions 

The first step in the analysis of our results was the obtaining of spatial profiles of 
several nebular parameters along the slit positions. The selected parameters were: c(H/3), 
ne, Te([N II]), re([0 III]), the intensity of several selected lines (H/3, C II 4267 A, O II 4649 
A, [O III] 4959 A, [Fe III] 4881 A, [N II] 5755 and 6584 A, [O I] 6300 A, and [S II] 6717 + 
6731 A), the 0^+/H+ ratio obtained from CELs and RLs, and the C"'""'"/H"'" ratio obtained 
from RLs. In figures 4 to 8 we show some selected spatial profiles of slit positions 3, 6, 1, 4, 
and 5, respectively. Slit positions 3 and 6 are the most interesting ones and we will focus our 
discussion on their main features. Slit position 3 crosses the Orion bar and the Herbig-Haro 
(H-H) objects HH 203 and HH 204 and slit position 6 passes through the brightest part of 
the nebula at the southwest of the Trapezium cluster, two proplyds: 159-350 and 177-341, 
and HH 530. 

The spatial profiles of show a large range of variation across the slits, with local 
maxima associated with the position of proplyds, H-H objects, the Orion bar, and the bright 
zone at the southwest of the Trapezium (see figures 4a to 8a). The highest densities are 
found at the proplyd 159-350 that has been observed in slit positions 5 and 6 (see figures 
8a and 5a, respectively). This object shows of the order of 6 x 10^ and 2 x 10^ cm~^ in 
positions 5 and 6, respectively; whereas the proplyd 158-326 -near 9^ Ori C- shows values 
somewhat larger than 4 x 10^ cm~^ (See Figure 6a). The densities at the brightest zone 
of the nebula -the southwest of the Trapezium- are about 2.5x10^ cm~^ (see Figure 5a). 
It is obvious that so high n^, determinations based on the [S II] doublet are not totally 
confident because they are at the high-density limit of this indicator. The H-H objects are 
also associated with local peaks of ne, but not as high as in the proplyds, in fact HH 202, 
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HH 203, and HH 204 show maxima between 6000 and 1 x 10^ cm~^. In Figure 4a, the 
dashed hne marks the obtained from the "whole sht" (integrated) spectrum -2460 cm~^- 
whereas the minimum and maximum values are about 700 and 8000 cm~^, covering a range 
of one order of magnitude. In Figure 5a, we can see that the range of variation of is also 
dramatic along slit position 6. In this figure we also compare the values corresponding to 
the "whole slit" spectrum -5500 cm~^- and the "background gas" -4700 cm~^- one, that 
corresponds to an integrated spectrum excluding the emission of the proplyds. As we can 
see, the measured density increases in 800 cm~^ -about a 17%- when we include the emission 
of the proplyds in the integrated spectrum of the slit. 

The spatial profiles of Te([N II]) and Te([0 III]) show very interesting features (see figures 
4b to 8b). The proplyds observed in positions 1, 5, and 6 show quite prominent spikes of 
Tc([N II]) and lesser or almost absent ones of Te([0 III]). In slit position 6 (Figure 5b) the 
Tc([N II]) increases locally about 70% at proplyd 177-341, a similar spike shows 170-337 in 
slit position 5 (Figure 8b). The Te([N II]) is also higher than the mean by 50% and 40% at 
the position of proplyd 159-250 in slit positions 5 and 6, respectively. In all the proplyds, 
the increase of Te([0 III]) is only of a few hundred K at most. In Figure 5b, we can also see a 
relatively broad -about 5" wide- spike of Te([N II]) at 130", where the ter aperature in c reases 
about 15%. This feature is not related to any local structure reported by lBally et al.l (120001 ) 
but with a conspicuous dark globule that can be seen at the edge of the bright zone at the 
southwest of the Trapezium. There are also less important temperature spikes related to 
some H-H objects. In Figure 4b, we can see that Te([N II]) increases about 15% at HH 204 
and that the increase is slightly lower in the case of Tc([0 III]). In contrast, HH 203 do not 
show that behavior and its temperatures are similar to those of the surrounding gas. In 
HH 202 -slit position 4- the increase of Tc([N II]) and Tc([0 III]) are only about 8% and 
5%, respectively (Figure 7b). The other H-H objects we have observed -HH 529 and HH 
530- do not show temperature variations with respect to the surrounding gas (figures 6b 
and 5b). The different behavior of the temperatures in the H-H obj ects does not seem t o 
be correlated with the velocity of their associated flows as reported by lHenney et al.l (120071 ). 
There is a last interesting feature regarding the temperature profiles that can be seen in 
Figure 4b. Although Te([N II]) is almost always a few hundred K larger t han Tp([0 III]) 



in all the slit positions (this has been also observed in previous works, e.g. [Baldwin et al. 



199ll : lRubin et al.ll2003h . the zone around the Orion bar -between 130" and 150" in Figure 
4b- shows a reversal of this relation just in the inner 10" of the bar. In contrast, Te([N II]) 
shows a local increase just outside the bar. This local increment of Te([0 III]) we see in this 
particular zone could be related with the highly ionized jet -specially bright in [O III]- that 
leads to HH 203 and HH204 (see lDoi et al.ll2004j ). In fact, the posit ion of this z one c oincides 
with a knot of high [O III] emission that can be seen in figure 8 of iDoi et al.l (120041 ) that is 



- 16 - 




Fig. 4. — Spatial profiles of several nebular parameters along slit position 3. Positional 
measurement along the slit goes from southeast to northwest (see Figured]). The position of 
the Orion bar and the Herbig-Haro objects HH 203 and HH 204 are indicated. Typical error 
bars are included in some of the diagrams, (a) ne, the horizontal long-dashed line gives the 
average value of for the "whole slit" (integrated) spectrum, (b) upper pannel: Te([0 III]) 
(dashed-dotted line) and Te([N II]) (continuous line); lower pannel: re([N Il])/Te([0 III]) 
ratio, the horizontal continuous and long-dashed lines represent the value of that ratio for 
the "background gas" and 'whole slit" spectra, respectively, (c) F([Fe III] 4881). (d) upper 
pannel: F{C II 4267) (continuous line) and F{0 II 4649) (dashed-dotted hne); lower pannel: 
F(C II 4267)/-F(0 II 4649) ratio, the vertical grey lines indicate zones without a reliable 
measurement of the O II 4649 line, (e) upper pannel: F{\0 III] 4959) (continuous line) 
and F(0 II 4649) (dashed-dotted line), lower pannel: F([0 III] 4959)/F(0 II 4649) ratio, 
the vertical grey bands indicate zones without a reliable measurement of the O II 4649 line, 
(f) upper pannel: 12 -|- log(0"'"^/H"'") determined from RLs (continuous line) and CELs 
(dashed-dotted line); lower pannel: ADF(0+"'"), the straight continuous and long-dashed 
lines represent the value of the ADF for the "background gas" and 'whole slit" spectra, 
respectively, the vertical grey bands indicate zones without a reliable determination of the 
0^^/H+ ratio from RLs. 
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Fig. 5. — Spatial profiles of several nebular parameters along slit position 6. Positional 
measurement along the slit goes from northeast to southwest (see Figured]). The position 
of the proplyds 177-341 and 159-350 and the Herbig-Haro object HH 530 are indicated. 
Typical error bars are included in some of the diagrams, (a) ne, the horizontal continuous 
and long-dashed lines represent the value of Ue for the "background gas" and 'whole slit" 
spectra, respectively, (b) upper pannel: Te([0 III]) (dashed-dotted line) and Te([N II]) 
(continuous line); lower pannel: Te([N Il])/Te([0 III]) ratio, the horizontal continuous and 
long-dashed lines represent the value of that ratio for the "background gas" and 'whole slit" 
spectra, respectively, (c) -F([[N II]] 5755) (continuous line) and -F([[N II]] 6584) (dashed- 
dotted line), (d) upper pannel: F{C II 4267) (continuous line) and F{0 II 4649) (dashed- 
dotted line); lower pannel: F(C II 4267) /F(0 II 4649) ratio, (e) upper pannel: F([0 III] 
4959) (continuous line) and F{0 II 4649) (dashed-dotted line); lower pannel: F{[0 III] 
4959)/F(0 II 4649) ratio, (f) upper pannel: 12 + log(0++/H+) determined from RLs 
(continuous line) and CELs (dashed-dotted line); lower pannel: ADF(0"'"^), the horizontal 
continuous and long-dashed lines represent the value of the ADF for the "background gas" 
and "whole slit" spectra, respectively. 




Fig. 6. — The same as Figure [5] but for slit position 1. Positional measurement along the 
slit goes from north to south (see Figured]). 
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Fig. 7. — The same as Figure H] but for slit position 4. Positional measurement along 
slit goes from southeast to northwest (see Figured]). 
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Fig. 8. — The same as Figure [5] but for slit position 5. Positional measurement along the 
slit goes from northeast to southwest (see Figured]). 
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crossed by our slit (their box with the position-based identifier 209-446). 



Rubin et aLl (120031 ) obtained HST/STIS spectroscopy of a sht position very similar to 
our position 3 -their slit 4- and show its Tc([0 III]) spatial profile in their figure 3a. If 
we compare that figure with our Figure 4b we can see that the point-to-point dispersion of 
the temperature is substantially lower in our d ata. I n fact, whereas Te([0 III]) varies from 



6500 to 12000 K in the slit 4 of iRubin et all (l2003f ) the variations are only from 8000 to 



9000 K in our slit position 3. This fact has two possible explanations: a) the presence of 
real temperature variations with typical spatial scale between O'.'S (the spatial resolution of 
Rubin et al. data) and 1"2 (our res o lution ) in the plane of the sky, or b) the temperature 
variations reported by iRubin et al. J2003h are spurious and produced by the much lower 



signal-to-noise ratio of their data. We think that the second explanation is perhaps the most 
likely one considering that the deepest exposures obtained by Rubin et al. for the spectral 
range containg [O III] 4363 A were about 1060 s and our exposures were between 1800 and 
2200 s long. Moreover, an important factor should be added due to the different apertures 
of the telescopes used in our and Rubin et al.'s datasets -4.2m and 2.4m, respectively- as 
well as the smallest element of spatial resolution used in Rubin et al.'s observations. 

It is clear that the behavior of Tc([N II]) and Te([0 III]) at the positions of prop 1yds and 
H-H objects is different. In figures 4b to 8b we also include the ratio of both temperatures 
showing clearly that the increase of Te([N II]) is higher than that of To([0 III]) in proplyds. 
The presen ce of Tp([N Ilj) enha ncements in the proplyds of the Orion nebula was previously 



reported by lRubin et al.l (120031 ) and they interpret this as the effect of coUisional deexcitation 
on the nebular lines of [N II] due to the high densities of these objects. We have explored 
that possibility comparing the spatial profile of the intensity of [N II] 5755 and 6584 A lines. 
Both lines come from upper levels with very different critical densities -7.9x10^ cm~^ in the 
case of [N II] 5755 A line and 5.8x10^ cm~^ in the case of [N II] 6584 A line. In figures 5c 
and 6c, we can see that the two proplyds with the largest electron densities (ne > 2-6x10^ 
cm~^): 159-350 (also observed in Figure 8c) and 158-323 show a spike in the brightness of 
the [N II] 5755 A line, whereas the [N II] 6584 A line do not show a so clear increase in its 
brightness with respect to the emission of the surroundings. However, the proplyds showing 
the lowest electron densities (rie < 1x10"^ cm~^) 177-341 and 170-337 (figures 5c and 8c) 
show similar localized enhancements of the intensity of both [N II] lines, indicating that 
coUisional deexcitation seems to be not so important in these two proplyds. 

In order to further explore if the different behavior of the auroral and nebular Te([N II]) 
lines in some proplyds is due to coUisional deexcitation, we have constructed Figure M In 
this figure, we show the theoretical curves of the Ue and Te pairs that reproduce the observed 
range of values of the [N II] 5755/6584 line ratio in the proplyds, as well as the line of the 
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lower limit of corresponding to the lowest value of the [S II] 6731/6717 ratio measured in 
these objects. The theoretical predictions have been constructed with e missivities calculated 
by the photoionization code PHOTO as described in IStasihskal (120051 ) . gently provided by 
Grazyna Stasihska. Unfortunately, the results of Figure are not conclusive, but considering 
that the estimated rze of the proplyds should be a lower limit of the true one, the permitted 
area of the diagram indicates that coUisional deexcitation should be acting. Moreover, if 
that is correct, the true electron temperature should be lower than that actually indicated 
by Te([N II]), but perhaps not too different that the values corresponding to the surrounding 
ionized gas (never lower than 3000 K). It is important to note that the re([N II]) values we 
obtain for the "whole slit" and the "background gas" spectra of positions 5 and 6 only show 
differences of the order of a few ten K. This indicates that the contribution of proplyds in 
the integrated spectrum does not produce a substantial increase of the electron temperature 
derived. On the other hand, in contrast with what happens in proplyds, the increase of 
Te([N II]) that we see in H-H objects does not seem to be related to collisional deexcitation. 
This is suggested because the H-H objects show lower values of rig and lower Te([N II]) peaks 
than the proplyds, their Te([N II]) and Te([0 III]) show similar enhancements (see figures 4b 
and 7b), and the spatial coincidence of conspicuous similar peaks in both [N II] 5755 and 
6584 A hues. Therefore, the [N II] emission of the H-H objects should not be affected by 
substantial collisional deexcitation and their electron temperature is higher probably because 
the action of an additional source of heating, perhaps related to shock excitation. 



5.2. Line Fluxes 

One of the main spectral properties of H-H objects is their strong emission in [Fe III] 
lines. In Figure 4c, we can see the spatial profile of F([Fe III] 4881)//(H/?) line ratio along 
slit position 3. In HH 204, the intensity of that line is a factor 5 brighter than in the 
Orion bar. That object also shows a similar enhancement in the [O I] 6300 A line and more 
moderate ones in the [S II] and [N II] lines. In the case of HH 202 -observed in slit position 
4-, the intensity of [Fe III] 4881 A line increases by a factor 10 (Figure 7c). The rest of the 
H-H observed: HH 204, HH 529, and HH 530 also show enhancenments in [Fe III], [O I], 
[N II], and [S II] lin es. These are common spectral features in shock-excited objects (see 
Hartigan et al. 1987 ). In contrast, the proplyds show [O ij^ spikes of moderate intensity and 



reversed spikes in the [S II] 6717, 6731 A lines, perhaps also due to collisional deexcitation 
because of the rather low critical densities of the upper levels of those [S II] transitions. The 



''The [O I] 6300 A line could not be corrected from sky emission contribution. In any case, spatial 
variations of the intensity of this line along the slit can only be due to intrinsic nebular variations. 
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Fig. 9.— Te vs. 



rip 



The continuous lines represent the theoretical curves of the rig and Te 



pairs that reproduce the observed range of values of the [N II] 5755/6584 line ratio measured 
in the proplyds. The dashed line represents the lower limit of the Ue corresponding to the 
lowest value of the [S II] 6731/6717 ratio measured in the proplyds. The grey area represents 
the permitted zone for the average spectral properties of the proplyds. The theoretical 
predictions have been c onstructed with emissivities calculated by the photoionization code 
PHOTO as described in IStasihskal (120051 ). gently provided by Grazyna Stasihska. 
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[O I] emission emerges from the photodissociation of OH in t he H/H2 front that hes close to 
the protoplanetary disk surface (IStorzer fc Hollenbachlll998l ). 

In figures 4d to 8d, we show the spatial profiles of the pure RLs of C II 4267 A and O II 
4649 A along the slit positions. The spatial distribution of the C II and O II lines are very 
similar. However, in the cases of the slits positions passing through the center -specially in 
positions 1 and 6- there is a slight decrease of the C II/O II ratio towards the central parts 
of the nebula. This variation could be due to the increase of the C^^ ionization fraction near 
the Trapezium stars. 

In figures 4e to 8e, we show the spatial profiles of the O II 4649 A and [O III] 4959 
A lines, that show a fairly similar spatial distr ibution in all the slit positions. This behav- 



ior is very different to t hat observed in PNe (iLiu et al.l l200d : iGarnett fc DinersteinI 12001 



Krabbe fc Copettil l2006l ) where the O II line emission peaks closer to the central star than 



the [O III] line. We only find some localized enhancements of the [O IIl]/0 II ratio, which 
are related to the position of some proplyds (see figures 5e and Be). These enhancements of 
the [O III] /On ratio are correlated with the increase of the continuum due to the emission 
of proplyds that show a large number of absorption features that produce a decrease in the 
intensity of the O II 4649 A line (see i CT) . 



5.3. The Abundance Discrepancy Factor 

Finally, figures 4f to 8f show the spatial variation of the 0"'"^/H"'" ratios obtained from 
CELs and RLs as well as the AD factor, (hereinafter ADF) defined as: 

ADF(0++) = log(0++/H+)RL, - log(0++/H+)cELs, (3) 

for all the slit positions. The most interesting result concerning these figures is that the 
ADF remains fairly constant along most of the zones of the nebula observed, showing values 
between 0.15 and 0.20 dex, in agreement with the determinations by Esteban et al. (1998 
and 2004) based on deep echelle spectrophotometry of selected small areas of the Orion 
nebula. It is necessary to note that the behavior of the ADF close to the proplyds is not 
confident because the O II lines are not properly measured in most of these objects due to 
the aforementioned strong increase of the continuum that makes difficult the measurement 
of weak lines and the effect of possible absortion features of their spectra. In any case. Figure 
5f shows a quite convincing decrease of about 0.10 dex in the 0^+/H+ ratios obtained from 
CELs at proplyd 177-341. 

Another remarkable feature of the ADF along the slits can also be seen in Figure 5f, 
where we find a slightly higher ADF in an area between proplyds 159-350 and 177-341. It 
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is interesting to note that slit position 5 (Figure 8f) also shows the same higher values of 
the ADF between proplyd 159-350 and the area just at the north of 177-341. That common 
behavior in both slit positions indicates that the local increase of the ADF should be real. 
This enigmatic zone of relatively high ADF has an apparent diameter of about 30" and is 
located about 23" south of 6^ Ori C star. There is not an apparent morphological and/or 
kinematical feature related to this zone. 

In Figures 5f to 8f we also include the values of the ADF corresponding to the "whole 
slit" and the "background gas" spectra. It can be seen that the ADF of the "background 
gas" is consistent with the average value of the different individual apertures extracted from 
the slit positions. However, the ADF obtained for the "whole slit" spectrum -which is almost 
coincident with that of the "background gas" in slit position 1, 3, and 4- is substantially lower 
-about 0.1 dex- in slit positions 6 and 5 (figures 5f and 8f, respectively). The reason of this 
surprising result is most probably related to the strong contribution of the proplyds to the 
continuum of the integrated spectra and that this contribution is producing some absorption 
in the O II lines. In fact, in slit positions 5 and 6, the continuum adjacent to the O II 4649 A 
line is about a factor 1.7 higher in the "whole slit" spectra than in the "background gas" one 
(see Figure [TU] for sht position 5). This higher continuum is also associated with a decrease 
of the intensity ratios of O II lines with respect to H/5. For example, the /(O II 4649)/J(H/5) 
ratio is about a 25% lower in the "whole slit" than in the "background gas" spectrum of slit 
position 6. 

In figures 4f and 7f, we can note that the ADF is particularly large -reaching even values 
up to 0.6 dex- along HH 203, HH 204 and HH 202. The large increase of the ADF is due to 
the low values of the O+^/H'^ ratios determined by CELs in these zones. In contrast, the O^"*" 
abundances determined by RLs do not show so strong localized decrease. It is interesting 
to compare the behavior of the ADF at the H-H objects and at the Orion bar in Figure |H 
In the bar, we can see a similar decrease of the 0++/II+ ratios determined from both kinds 
of lines but producing an ADF similar to the mean value along the slit. As it has been 
discussed above, coUisional deexcitation does not seem to affect the intensity of the nebular 
[N II] lines at HH 203 and HH 204 and, therefore, considering the larger critical densities 
of the [O III] lines, this effect is even less likely to be producing the observed decrease of 
the 0^"'"/H"'" ratios determined from CELs. The observed behavior of the ADF in the H-H 
objects could be explained because the presence of a localized heating due to a non-radiative 
process -most likely shock excitation- that would lead to the derivation of a lower 0++/H+ 
(CELs) ratio, and this seems to be the case for HH 204 that shows a conspicuous spike of 
Te([0 III]) (Figure 4b). However, this hypothesis fails to reproduce the presence of high 
ADFs in HH 203 and HH 202, where there are not spikes of Te([0 III]). Another explanation 
is the presence of a localized high in those particular zones but, unfortunately, we can not 
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check this possibihty with our data. Finally, the other H-H objects observed -HH 529 and 
HH 530- do not show any distinguishable enhancement of the ADF. 



6. Radial Distributions of Some Relevant Nebular Parameters 

We have constructed radial distribution diagrams of several nebular parameters combin- 
ing the data of the different slit positions and projecting the angular distance of the center 
of each individual ID spectrum with respect to 6^ Ori C (Figure [TTj) . In these diagrams, 
we have only included the data of those apertures corresponding to the "background gas", 
i.e. we have excluded all the points associated with proplyds or H-H objects. We include 
the same points in the three diagrams of Figure [11] and only those for which the ADF was 
calculated. These diagrams permit to study possible large scale variations of the properties 
of the ionized gas across the nebula. 

In Figure 11a, we show the radial distribution of Te([0 III]), that shows a slight but 
clear general decrease with increasing distance from 6^ Ori C. In this figure, it can be seen 
a bump in the temperature distribution of slit position 3 between 90" and 120". This bump 
corresponds to the local increase of Te([0 III]) that occurs at the inner part of the Orion bar 
(see also Figure 4b). We have made a least-squared linear fit to the data shown in Figure 
11a but excluding those points belonging to the Orion bar. The result is: 

Te([OIII])(K) = 8540 -3.6 X r, (4) 

where r is the distance from 6^ Ori C in arcseconds. The uncertainty of the slope is 0.2 
K arcsec"^ and the Spearman's correlat ion coefficient is 0.6 3 (hereinafter all the correlation 



coefficents we use are Spearman's ones). IWalter et al.l (119921 ) obtained a radial distribution of 



Te([0 III]) with a positive slope of 6.7 K arcsec"^ in the northwestern quadrant of the nebula, 
fitting their own data and other from the literature. The result of these authors disagrees 
with the clear behavior shown by our data in Figure 11a, that seems to be independent of the 
location of each slit position an d even includes zon es of the northwestern quadrant -our slit 



position 4. On the other hand, lO'Dell et al.l (120031 ) do not find significant spatial variations 



of T'e([0 III]) across the nebula in their high- resolution map of the [O III] ratio obtained 
with the WFPC at the HST. As we can see, the different results about the large-scale 
temperature variations in the Orion nebula are contradictory, however, we consider that our 
dataset is more reliable than the previous ones because it is based on homogeneous higher 
signal-to-noise ratio spectroscopic observations. 



In Figure lib, it is evident that Te([N II]) also shows a radial decrease. The least-squared 
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Fig. 10. — Comparison of the spectral area around the multiplet 1 of O II for the "whole 
slit" (up) and "background gas" (down) spectra of sht position 5. Note the rather different 
contribution of the continuum emission. The "whole slit" spectrum includes the emission of 
proplyds and the "background gas" spectra do not. 
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Fig. 11. — Radial distribution of the values of several nebular parameters along the Orion 
nebula. We have included the data of the apertures of all slit positions excluding those 
associated with proplyds or Herbig-Haro objects, i.e. the figures show the behavior of the 
featureless ionized gas. The abscissa axis indicates the projected angular distance of the 
center of each individual ID spectrum with respect to 9^ Ori C. Symbols are described in 
(a). Typical error bars are included in the different diagrams, (a) Te{[0 III]), the continuous 
line represents a least-squared linear fit to the data, (b) Tc([N II]), the continuous line 
represents a least-squared linear fit to the data, (c) ADF(0"'""'"), the dashed line represents 
a least-squared linear fit to the data, the continuous line represents a divided least-squared 
linear fit to the data into two regions inside and outside 40". The fits to the data are given 
in the text. 



- 29 - 



linear fit to the data gives: 

Te([N II]) (K) = 9460 - 7.8 xr . (5) 
We can see that the slope of the fit is larger than that obtained for Te([0 Ill j) and, in this 



case, t he result agrees qualitatively with the temperature gradient obtained by IWalter et al. 



(I1992I ) from the [N II] lines (—17.1 K arcsec"^) in the inner 170" of the nebula. The un- 
certainty of the slop e of our fit is 0.6 K arcsec"^ and its correlation coefficient is 0.52. 
Sanchez et al.l (120071 ) obtain a bidimensional map of the spatial distribution of Te([N II]) 
finding also that this parameter is higher near the Trapezium stars and drops towards the 
outer zones of the nebula, in qualitative agreement with our result. 

Finally, in Figure 11c, we show the radial distribution of the ADF(0"'"+) across the 
nebula. In this figure, part of the relatively large dispersion is due to the large observation 
uncertainty in several apertures. In fact, most of the points with the highest and lowest 
values of the ADF correspond to low signal-to-noise ratio determinations at the edges of slit 
positions 1, 4, and 5. Despite the dispersion is relatively large, the first visual impression 
suggests that the ADF is rather constant or very slightly decreasing towards the outer parts 
of the nebula. In fact, the least-squares linear fit of all the points gives: 

ADF(0++) = 0.2050 - 0.0005 x r , (6) 

but with a very low correlation coefficient (0.32). However, other possibilities of fitting are 
possible. A more detailed inspection of Figure 11c shows that the ADF seems to increase 
toward the center in the innermost 40" of the nebula and that this parameter becames 
basically constant for larger distances (the inner points of slit positions 1 and 5 consistently 
suggest the same tendency). We have divided the least-squared linear fit to the ADF into 
two regions, the first for the points in the inner 40" and the second for the points outside 
this area -for which we have assumed a zero slope. The fits are the following: 

ADF(0++) = 0.279 - 0.003 x r (r < 40") , (7) 

with a correlation coefficient of 0.58, and 

ADF(0++) = 0.163 (r > 40") . (8) 

This apparent increase of the ADF toward the inner zones of the nebula can be ultimately 
related to two other tendencies we have found in our data. Firstly, the clear correlation 
illustrated in Figure 11a showing that To([0 III]) is higher in the zones near the Trapezium 
stars and, secondly, the possible weak correlation between the ADF and Te([0 III]) that will 
be discussed in J|71 Therefore, the ADF seems to increase very slightly in the inner and 
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systematically hotter zones of the nebula. If this behavior is real, it would indicate that 
whatever process is producing the ADF it increases somehow with the local ionization of the 
gas or that it depends on the distance between 6^ Ori C and the ionization front, from where 
most of the nebular emission comes. This last possibility arises because the blister geometry 
of the Orion nebula, the smaller distance of the Trape zium stars to the ma. in ionization front 
occurs precisely behind that massive star cluster (see lWen fc 0'Delllll995l ). 



7. Correlations Between the ADF and Other Nebular Properties 

In order to shed some light on the physical nature of the ADF problem, we have explored 
the relationship of this parameter with other nebular properties determined from our ID 
spectra. Figure [T2] illustrates the dependence of the ADF(0+"'") with respect to a) c(H/5), b) 
ne, c) Te([0 III]), d) Te([N Il])/Te([0 III]) ratio, e) 0++ abundance determined from CELs, 
and f) 0"'"+/0"'" ratio. As in Figure [TT| we have only included the data of the apertures 
corresponding to the "background gas" . Figure 12a indicates that the ADF does not depend 
on the amount of dust present in the line of sight of each individual ID spectrum. Also, 
there is no correlation between the ADF and ne as it shown in Figure 12b. However, from 
this figure it would seem that ADFs larger than about 0.15 dex are only found in zones of 
ric < 1.3 X 10^ cm~^. This apparent trend is incidentally produced by the apertures of slit 
position 1, which are located near 6^ Ori C and show the largest values of ADF, a behavior 
that does not seem to be related to Ue but perhaps with the increase of Te([0 III]), as we 
will see below . The large dispersion of the ADF for the lowest densities shown in Figure 
12b is an observational bias. The less dense zones show the lowest surface brightnesses and 
therefore the faintest spectra, implying larger uncertainties in their ADF determinations. 
The behavior of the ADF with respect to ^^^([O III]) is shown in Figure 12c and seems 
to show an apparent weak positive correlation, although Te([0 III]) shows a rather narrow 
interval of values in the nebula. 

The other temperature indicator, Te([N II]), shows a similar behavior and this is reflected 
in Figure 12d, where we represent the ADF versus the Te([N II])/Tc([0 III]) ratio. It is 
interesting to note that in Figure 12d the data points are distributed delineating the cloud of 
uncertainty around the average value of the ADF and Te([N II])/Te([0 III]), illustrating that 
the portion of the nebula observed shows a tight proportionality between both quantities. 
The lack of correlation between the ADF and the ratio of electron temperatures of ions 
located in different zones of the nebula indicates that the natural temperature gradients 
that should exist in ionized nebulae -due to the different spatial location of the main nebular 
coolants- do not play a significant role in producing the observed ADF in the Orion nebula. 
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Fig. 12. — The abundance discrepancy factor of ADF(0+''"), vs. several nebular 

parameters. We include the data from all sht positions. Symbols are described in (b). 
Typical error bars are included in the different diagrams, (a) c(H/3). (b) Ug. (c) Te{[0 III]), 
(d) re([N II])/Te([0 III]) ratio, (e) 0++ abundance determined from CELs. (f) 0++/0+ 
ratio. 
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This results is consistent with that found by lGarcia-Rojas fc EstebanI (120071 ) in their analysis 
of integrated spectra of a sample of Galactic and extragalactic H II regions. Figure 12e also 
shows a lack of correlation between the ADF and the 0~^~^ /'H~^ ratio obtained from CELs 
-a similar diagram is found in the case of the 0"*"+ abundance obtained from RLs, although 
the range of variation of the 0"'""''/H"'' ratio is very narrow. Only the points belonging to 
the Orion bar show values of the O^'^ abundance between 0.2 and 0.3 dex lower than the 
typical ones of the rest of the nebula, but their corresponding ADF is similar. Finally, Figure 
12f represents the behavior of the ADF with respect to the 0++/0+ ratio and it is evident 
that, despite the relatively large range of values covered by the ionic ratio, there is not an 
apparent trend between both quantities. 



8. Temperature Fluctuations 



Following the formulation proposed by iPeimbertI (Il967l ). the temperature fluctuation 
over the observed volume of a nebula can be parametrized in terms of the average tempera- 
ture. To, and the mean-square electron temperature fluctuation, t^, defined as: 



To 



/ neUidV 



and 



/ (Te - ToYn^mdV 



(9) 



(10) 



T2 / n,n,dV 

where nj is the ion density. The integrations are calculated over the entire volume, and the 
element of volume, dV, can be expressed as dldA, the product of the elements of length of 
column along the line of sight and surface area in the plane of the sky, respectively. Our 
spatially resolved spectroscopical data do not permit to obtain a direct determination of 
along the line of sight. However, a discrete estimation of in the plane of the sky -t\- can 
be obtained through the point-to-point determinations of the Tg we have obtained from the 
i ndivi d ual ID spectr a extra cted a long the slit positions. Fo llowing a similar procedure as [Liu 
Jl998h . iRubin et al.l J2003h . and iKrabbe fc Copettil J2005h and assuming that Ue ^ n(H+) 
in all the points of the nebula, we can compute Tq a and t\ using the following equations: 



E,Te,ng,[n(X+0/n(H+)], 

E,<,h(x+0/^(H+)], 



and 



E,(7;.-To.a)X.[^(x^')/^(h+)L 



(11) 



(12) 
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where T^j, rigj, and [ra(X+')/r2(H+)]j are the electron temperature, electron density, and the 
ionic abundance of the species in the j-th aperture extracted from a given slit position. 
The Tej used in the equations above corresponds to average temperatures along the line of 
sight that crosses the nebula at each j-th aperture and it can be expressed as: 



T ■ 



(13) 



Taking into account this fact, it is likely that the t\ we obtain is substantially lower than 
and, more stric tly speaking, it should be considered a lower limit to (see further argumen- 
tation given by lRubin et al.ll2003l : lO'Dell et al.ll2003l ). In Table HI we summarize the values 



of t\ we obtain for each slit position and for the two ions for which we have determinations 
of Te. We have included all the points with Te determination available in the sums, even 
those belonging to proplyds or H-H objects in order to explore the effect of these structures 
into the derivation of t^. Part of the value of t\ we calculate for each ion comes from errors 
in the measurement of the emission line ratios, so the intrinsic t\ must be corrected by the 
relative r nean quadratic error of the Tp measuremen ts, t\^^, with the simple relation t\ — 
t\^^ (see lO'Dell et al.l l2003l : iKrabbe fc Copettil 120051 ) . The values of t\^^ for each ion and 



slit position are also included in Table H] and are always lower than the corresponding t\ 
except in the case of slit position 3, where the errors are slightly higher, indicating that 
most of the temperature variation along slit position 3 is produced by measurement uncer- 
tainties. The corrected t\ we obtain is very low in all the cases, t^fO" ^"*") ranges from ~0 
to 0.0011 and ti(N+) from ~0 to 0.0112. However, Ip.ubin et al.l (EooJ obtain t\{0++) = 
0.0068-0.0176 and ti(N+) = 0.0058-0.0175 for the different slit positions they observe in 
the Orion nebula, values which are relatively consistent with ou rs in the case of t^4 (N "'") but 
considerably larger in the case of t\{0~^~^). On the other hand, lO'Dell et al.l (l2003n obtain 
^ilO^^) = 0.0050-0.0156, values in good agreement with those obtained by iRubin et al. 



(120031 ) but also much higher than our determinations. The reason of this discrepancy is 
difficult to ascertain. The resolution element of each set of observations is rather different, 



our apertures are of 1'.'2 x 1'.'03, and those of iRubin et all (120031 ) and lO'Dell et all (120031 ) 
are 0'.'5 x 0'.'5 and O'.'l x O'.'l, respectively. If the spatial resolution has something to do 
with the different t\{0~^^) obtained by us and the previous works, this would imply that the 
small-scale spatial variations of Te([0 III]) have a characteristic size between 0'.'5 and I'.'O and 
that those related to Te([N II]) are larger, at least of the order of the size of our apertures. In 
any case, this seems to us a rather unlikel y sce nario. Anoth e r diff erence between our deter- 
minations and those of iRubin et al.l (120031 ) and lO'Dell et al.l (120031 ) is that those authors use 
a relation for determining from the [O III] line ratio that is valid in the low-density limit, 
while our calculations are based on the solving of the statistical equilibrium equations for 
the He measured for each particular aperture. We have explored the effect of this different 
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procedure applying equation 5 of iRubin et al.l (120031 ) -also used by lO'Dell et al.l (120031 )- to 
our data of slit position 6. In Figure [13], we show the ratio of the Te([0 III]) computed 
using the low-density limit approximation equation and our own determinations finding a 
systematical bias of about 1.12 and a small bump (~5%) coincident with the zone of the 
highest density at the southwest of the Trapezium (see Figure 5a). The use of both different 
procedures do not substantially increase the Tg dispersion. In fact, using the low-density 
limit approximation, we obtain t^(0"^^) = 0.0008 for slit position 6, only 60% larger than 
our calculations. Therefore, the use of a low-density limit approximation does not explain 
the large differences of t\{0^^) find by us and the other authors. A last possibility is the 
different signal-to-noise ratio of the datasets we are comparing. As it was commented in 



^ 15.11 this is a likely source of discrepancy in the case of the long-slit data of iRubiri et al. 



(120031 ) ■ but this is more difficult to ascertain in the case of the data of lO 'Dell et al.l (120031 ). 
These authors estimate a representative probable error of about 4.2% for their point-to-point 
Te([0 III]) determinations, which corresponds to a t^g^ of about 0.0017, larger than our ob- 
servational errors, but much lower t han their norn i nal v alues of t\{0~^^), indicating that the 
temperature variations obtained by lO'Dell et al.l (120031 ) should be real if the errors are not 
largelly underestimated. Only observations combining very high signal-to-noise ratio and 
very high spatial resolution would ascertain this puzzle. 

As it has been stated before, our spatially-resolved spectroscopical data do not permit to 
determine along the line of sight -that we will denote as tf- but an indirect estimate can be 
obtained assuming that the ADF is produced by the presence of such temperature variations. 
Although this is still a controversial possibility, there are some pie ces of evidence indicating 
that t his may be co r rect, at least in the case of H II regions (see iGarcia-Rojas fc Esteban 
20071 ). lO'Dell et al.l (120031 ) show that the relation between t\, tf, and the total t"^ in three 
dimensions is: 



tl+ <tf> 



"A^ - - , (14) 

where <tf > is the average over all lines of sight. In Table HI we include the <t^(0^+)> we 
obtain for each slit position. These values have been estimated from the average ADF(0^+) 
of the individual apertures of each slit position. From the table, it can be seen that the 
<t^(0"'"^)> is rather similar for the diffe rent slit positions and co nsistent with previous 
determinations -between 0.018 and 0.028 (lEsteban et al.lll998l . |2004| )- except in the case of 
slit position 3, that shows a rather higher value. This large temperature fluctuation comes 
from the also larger ADFs we obtain at HH 203 and HH 204 (see Figure 4f). Taking into 
account that the values of t^(0"'""'') are much lower than < t'f(0~^~^)>, it c an be considered 



that t2(0++) ^ <tf (0++)>. Following the argumentation of lO'Dell et al.l (l2003f ). this result 



would indicate that the hypothetical thermal inhomogeneities producing should be small- 
scale ones and unresolved by our data, i.e. smaller than our spatial resolution limit of about 
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Table 3. Selected lines 



A (A) 


Ion 


Mult. 


A (A) 


Ion 


Mult. 


4zD7.io 


TT 


6 


5517.71 


[r^^ TTTl 
[Ol lllj 




4340.47 


HI 


H7 


5537.88 


[CI III] 


IF 


4363.21 


[OIII] 


2F 


5754.64 


[Nil] 


3F 


4638.86 


on 


1 


6300.30 


[0 I] 


IF 


4640.64 


NIII 


2 


6312.10 


[S III] 


3F 


4641.81 


on 


1 


6363.78 


[01] 


IF 


4643.06 


Nil 


5 


6548.03 


[Nil] 


IF 


4649.13 


on 


1 


6562.82 


HI 


Ha 


4650.84 


on 


1 


6583.41 


[Nil] 


IF 


4661.63 


on 


1 


6716.47 


[SII] 


2F 


4711.37 


[Ar IV] 


IF 


6730.85 


[S II] 


2F 


4861.33 


HI 


H/3 


7135.78 


[Ar III] 


IF 


4881.00 


[Fe III] 


2F 


7319.19 


[0 II] 


2F 


4958.91 


[OIII] 


IF 


7330.20 


[Oil] 


2F 


5006.94 


[OIII] 


IF 


7751.10 


[Ar III] 


2F 



Table 4. Values of t\, t\^^, and tf 



Slit 



Position 


t\{0++) 


tL.(o++) 


ti(N+) 




<t2(0++)> 


1 


0.0013 


0.0002 


0.0127 


0.0015 


0.028 


3 


0.0004 


0.0006 


0.0010 


0.0011 


0.038 


4 


0.0008 


0.0003 


0.0055 


0.0016 


0.029 


5 


0.0004 


0.0002 


0.0126 


0.0019 


0.025 


6 


0.0005 


0.0002 


0.0024 


0.0017 


0.025 
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Fig. 13. — Ratio of th e Tpi\0 III]) c ompu ted using the low-density hmit (l.d.l.) approxima- 
tion equation used by iRubin et ahl (120031 ) and our own determinations along slit position 6. 
Positional measurement along the slit goes from northeast to southwest (see Figure [1]). 
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1". On the other hand, ou r spatial resolution is unable to resolve the size of 10 —10 cm 
that IStasihska et al.l (120071 ) derive for the metal-rich droplets they claim to be the most likely 
explanation of the AD. 



9. Conclusions 

We have studied the spatial distribution of a large number of nebular quantities along 
five slit positions covering different morphological zones of the Orion nebula. The resolution 
element of the observations was 1'.'2 x 1'.'03. The studied quantities were c(H/?), ne, Te([N II]), 
re([0 III]), the intensity of several selected lines (H/3, C II 4267 A, O II 4649 A, [O III] 4959 A, 
[Fe III] 4881 A, [N II] 5755 and 6584 A, [O I] 6300 A, and [S II] 6717 + 6731 A), the 0++/H+ 
ratio obtained from collisionally excited lines (CELs) and recombination lines (RLs), and the 
C^^/H"*" ratio obtained from RLs. The total number of apertures or ID spectra extracted 
was 730. We have been able to determine the 0+~*"/H"'" ratio from the faint RLs of this ion 
in a 92% of the apertures. 

The spatial distribution of shows a large range of variation -larger than an order of 
magnitude- across the nebula, with local maxima associated with the position of protoplan- 
etary disks (proplyds), Herbig-Haro objects, the Orion bar, and the brightest area of the 
nebula at the southwest of the Trapezium. The proplyds show quite prominent spikes of 
rc([N II]) and much lesser ones of Te([0 III]). This fact could be due to coUisional deexcita- 
tion on the nebular lines of [N II] because of the high densities of these objects. Herbig-Haro 
objects also show somewhat higher values of Te([N II]) but, in this case, the origin could 
be related to extra heating of the gas due to shock excitation. The spatial distribution of 
the O II 4949 A and [O III] 4959 A lines is fairly similar along all the slit positions, a very 
different behavior to that observed in planetary nebulae. We have found that the abundance 
discrepancy factor (ADF) of O"^"*" -the difference between the O"*""*" abundance determined 
from RLs and CELs- remains, in general, rather constant along most of the observed areas 
of the nebula, showing values between 0.15 and 0.20 dex. However, there are some localized 
enhancements of the ADF, specially at the position of the Herbig-Haro objects HH 202, HH 
203, and HH 204. 

The combined data of all slit positions indicate a clear decrease of Te([N II]) and 
rc([0 III]) with increasing distance from the main ionizing source of the nebula, 9^ Ori 
C. On the other hand, the radial distribution of the ADF shows a rather constant value 
across the nebula except at the inner 40", where the ADF seems to increase very slightly 
toward 9^ Ori C. 
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We have explored possible correlations between the ADF of and other nebular 
quantities, finding a possible very weak increase of the ADF for higher electron temperatures. 
There are not apparent trends between the ADF and c(H/?), Ue, the Te([N Il])/Te([0 III]) 
ratio, 0++ abundance, and the 0++/0+ ratio. 

Our spatially resolved spectroscopy allows to estimate the value of the mean-square 
electron temperature fluctuation in the plane of the sky, a lower limit to the traditional 
parameter. We find very low values in all cases, result that is in contradiction with 
previous estimates from the literature. Our results indicate that the hypothetical thermal 
inhomogeneities -if they exist- should be lower than our spatial resolution limit of about 1". 

It is clear that further studies on the T^, chemical abundances, and ADF distributions at 
sub-arcsec spatial scales are necessary in trying to disentangle (a) whether small spatial scale 
temperature fiuctuations and/or metal-rich droplets are really present in the Orion nebula 
and HII regions in general, and (b) the origin of AD problem and its possible relation with 
t'^ and other nebular properties. The observations needed for this task are very difficult even 
for ground-based large-aperture telescopes and, by now, unfeasible with the current space 
telescopes and their available instrumentation. 

We thank G. Stasihska and M. Rodriguez for their fruitful comments and help. We are 
grateful to the referee, Y. Tsamis, for his careful reading of the paper and his comments. 
This work has been funded by the Spanish Ministerio de Ciencia y Tecnologa (MCyT) under 
project AYA2004-07466. 
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